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摘要 本文 在 Zhang-Yu 工 作 的 基础 上 , Xf Caladine-English 折 板 吸 能 器 的 能 量 吸收 能 力 进 行 了 细致 分 析 . 我 们 利 
Maple 软 件 对 Caladine-English 折 板 吸 能 器 的 Zhang-Yu 吸 能 器 方程 编写 了 通用 程序 , 并 首次 求 得 了 吸 能 器 的 广 
SUM CERO 随时 间 的 变化 规律 . 在 此 基础 上 , 对 于 不 同 参数 变化 对 Caladine-English 折 板 吸 能 器 的 能 量 吸收 能 
力 的 影响 进行 了 比较 研究 , 其 中 为 了 刻画 结构 的 吸 能 能 力 , 定义 了 吸 能 因子 . 为 了 更 好 的 了 解吸 能 器 的 运动 受 参 
数 变化 的 影响 , 我 们 计算 了 吸 能 器 动力 学 的 相 图 . 经 过 大 量 的 数值 模拟 , 发 现 从 吸 能 能 力 上 看 , 4 折 板 吸 能 器 应 当 
是 “质量 敏感 "型 吸 能 器 ,完全 文 持 Zhang-Yu 的 结论 : 即 Caladine-English 折 板 吸 能 器 是 “质量 敏感 "型 吸 能 器 . 


关键 词 质量 敏感 性 , 吸 能 结构 , 初始 缺陷 , UERGER, Maple 
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Energy absorbing capacity analysis of the Caladine-English crooked plates energy absorber 


Sun Bohua*'Chen Pinyuan* 


*( School of Civil Engineering, Xi'an University of Architecture and Technology, Xi'an 710055, China) 
*( Institute of Mechanics and Technology, Xi'an University of Architecture and Technology, Xi'an 710055, China) 
*( Faculty of Mechanical Engineering & Mechanics, Ningbo University, Ningbo 315211, China) 


Abstract Based on Zhang-Yu's work, this paper analyzes the energy absorption capacity of the Caladine-English crooked 
plates energy absorber in detail. We used Maple software to write a general program for the Zhang-Yu energy absorber 
equation of the Caladine-English crooked plates energy absorber, and for the first time obtained the generalized coordinate 
(rotation angle) of the energy absorber over time. On this basis, a comparative study was made on the influence of 
different parameter changes on the energy absorption capacity of Caladine-English crooked plates energy absorber. In 
order to define the energy absorption factor, the energy absorption capacity of the structure was described. In order 
to better understand the effect of the movement of the energy absorber by parameter changes, we calculated the phase 
diagram of the energy absorber dynamics. After a lot of numerical simulations, it is found that the 4 crooked plates energy 
absorber should be a “mass sensitive" type of energy absorber in terms of energy absorption capacity, which fully supports 


Zhang-Yu's conclusion: the Caladine-English crooked plates energy absorber is “Mass sensitive" type energy absorber. 


Key words mass sensitivity, energy absorbing structure, initial imperfection, four-hinge model, Maple 


引言 力学 现象 . 结构 在 服役 状态 难免 会 受到 冲击 荷载 的 


冲击 碰撞 是 日 常生 活 和 工程 实际 中 广泛 过 到 的 。” 作用 , 为 了 避免 结构 受到 破坏 , 就 必须 考虑 冲击 能 量 
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的 吸收 , 由 于 金属 在 不 可 逆 的 塑性 应 变 阶 段 能 够 吸 
收 大 量 能 量 , 而 被 广泛 加 工 为 承受 碰撞 和 冲击 的 吸 
能 结构 或 者 吸 能 器 . 在 通常 的 弹 塑性 力学 中 , 由 于 讨 
论 的 都 是 准 静 态 问题 , 即 这 时 荷载 是 缓慢 地 施加 于 
材料 或 结构 上 , 相应 的 变形 响应 也 比较 缓慢 , 所 以 惯 
性 力 可 以 忽略 不 计 . 但 是 , 在 工程 实际 中 过 到 的 动态 
问题 , 在 突然 的 强 荷 载 作 用 下 , 材料 或 结构 的 变形 速 
度 很 快 , 这 时 就 材料 的 应 变 率 和 结构 惯性 的 影响 比 
较 重要 . 就 是 说 吸 能 器 或 者 结构 的 设计 需要 考虑 应 
变 率 和 惯性 对 于 吸 能 能 力 的 影响 上 29. 
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图 1 (a) 第 I 类 典型 结构 模型 (b) 第 [I 类 典型 结构 模型 


Fig. 1 (a)Type I typical structure model (b)Type Itypical structure 


model 


对 于 (图 DD) 中 的 三 种 典型 吸 能 结构 ， Calla- 
dine 和 Englishi1 指 出 材料 应 变 率 敏感 性 和 惯性 效应 
是 第 开 类 结构 速度 敏感 性 的 主要 原因 . 显然 , 深入 研 
究 第 工 类 结构 的 敏感 性 , 对 于 设计 能 量 吸收 结构 以 
及 确定 动态 模型 试验 的 尺度 率 十 分 重要 . 

Lee 和 Symonds 首 次 采用 刚 塑 性 模型 研究 了 梁 在 
强 动 载荷 下 的 动力 响应 , 并 提出 了 移行 鲍 的 概念 9 
此 后 , 围绕 梁 在 冲击 荷载 下 的 动态 响应 , 结构 冲击 动 
力学 有 了 长 足 发 展 P1-30. Booth 等 外 在 做 薄 钢 板结 构 
的 动力 试验 时 发 现 , 相 较 于 缩 比 尺寸 试 件 , 原型 结构 
的 变形 和 碰撞 时 间 更 大 . Calladine 等 (3 在 试图 解释 
上 述 现象 时 发 现 , 依照 准 静态 载荷 -位 移 曲 线 的 整体 
形状 , 能 量 吸 收 结构 大 致 可 以 分 为 两 种 类 型 : 第 I 类 
具有 一 条 相对 “平坦 ”的 载荷 - 位 移 曲线 , 而 第 开 类 结 
构 的 曲线 “快速 上 升 ” 随后 “急剧 下 降 ” 对 两 类 结构 
他 们 设计 了 两 组 试验 , 发 现 第 卫 类 结构 的 变形 对 碰 
撞 速 度 要 比 第 I 类 更 敏感 . 换 而 言 之 , 在 总 的 输入 动 
能 恒定 不 变 的 情况 下 , 由 高 速度 碰撞 得 到 的 最 终 变 
形 比 低速 度 碰撞 的 最 终 变 形 要 小 , 这 种 速度 敏感 现 
象 第 工 类 结构 比 第 I 类 结构 更 为 显著 . 

Yu 等 C5 加 根 据 两 个 物体 间 非 弹性 碰撞 的 经 典 
理论 提出 四 贸 链 模型 ,研究 了 落 锤 碰撞 速度 和 质 


量 比 与 碰撞 时 瞬时 能 量 损 失 之 间 的 关系 , 深刻 地 
揭示 了 “速度 敏感 性 ”的 实质 是 “质量 敏感 性 ” 同时 
强调 了 “初始 几何 缺陷 ?对 “速度 敏感 性 "的 显著 影响 . 
Tam 和 Calladine?7] 进 行 了 大 量 试验 , 发 现 折 板 响应 由 
两 相 组 成 . 第 一 相 是 折 板 轴 向 压缩 , 持续 时 间 短 暂 ; 
第 二 相 是 四 贸 链 模型 动力 响应 阶段 , 即 刚 杆 绕 塑性 
贸 转 动 阶段 . 下 文 将 会 沿用 这 一 表述 . 

Su 等 R83 在 Zhang-Yu 理 论 模 型 的 基础 上 考虑 了 
弹性 变形 ， 并 成 功 预测 到 了 峰值 荷载 . Karagiozo- 
va 等 B030 利 用 板结 构 动态 弹 塑性 屈曲 的 模型 , 综合 
虑 了 惯性 效应 和 材料 应 变 率 效 应 , 对 折 板 动力 行为 进 
行 了 细致 研究 . 他 们 理论 预测 的 结果 与 Tam 等 R72 试 
验 结果 大 致 相同 . Honig 等 63 基 于 弹 - 黏 塑性 模型 , 采 
用 有 限 元 软件 ABAQUS, 对 碰撞 作用 下 具有 初始 缺陷 
的 折 板 进行 了 模拟 . 模拟 结果 与 Tam 等 R71 试验 结果 
大 体 一 致 . 户 文浩 等 采用 有 限 元 软件 ABAQUS, i 
拟 了 两 类 典型 的 能 量 吸收 结构 在 铁 锤 动态 冲击 下 的 
力学 行为 , 他 们 发 现 : 开 类 结构 惯性 敏感 性 较 工 类 结 
构 显 著 ; 能 量 吸 收 结构 的 反 力 的 峰值 大 小 与 冲击 速 
度 相关 , 冲击 速度 越 大 , 结构 反 力 的 峰值 也 越 大 . 

从 目前 所 能 看 到 的 文献 , 对 于 第 HI 类 典型 结构 模 
型 的 分 析 工 作 , 以 Zhang-YuD 的 工作 具有 代表 性 . 由 
于 这 项 杰出 的 工作 发 表 于 30 多 年 前 的 1989 年 , 可 能 
限于 当时 的 数值 计算 能 力 , Zhang-Yu 外 没有 给 出 基本 
变量 转动 角度 90 的 时 间 变 化 规律 和 数值 曲线 . 由 
第 I 类 典型 结构 是 很 好 的 吸 能 结构 , 为 了 可 以 对 
种 吸 能 结构 的 吸 能 机 理 进 一 步 的 了 解 , 所 以 有 必要 
在 Zhang-YuD 工 作 的 基础 上 , 利用 现代 的 计算 手段 对 
其 进行 更 加 细致 地 分 析 和 解剖 特别 是 希望 可 以 
对 Zhang-YuD 的 一 个 观点 , 即 认为 第 I 类 典型 结构 不 
是 “速度 敏感 性 ”的 , 而 是 “质量 敏感 性 ”的 科学 判断 ， 
做 出 定性 到 定量 的 科学 支持 . 

本 文 在 Zhang-YuD 工 作 的 基础 上 , 运用 工程 计算 
软件 Maple, 首次 求解 了 非 线 性 的 Zhang-Yu 吸 能 器 技 
术 方 程 , 得 出 了 第 工 类 结构 在 落 锤 动 态 冲 击 下 的 不 
同时 刻 的 折 角 与 角速度 的 动力 学 . 通过 与 Tam 等 中 
的 试验 结果 比较 ， 验 证 了 Zhang-YuD 中 模型 的 合理 
PE, 同时 总 结 出 折 板 动力 响应 第 二 相 一 一 刚 杆 绕 贸 
转动 阶段 的 一 些 规 律 . 具体 讲 , 在 第 2 节 , 给 出 了 求 
解 Zhang-Yu 方 程 的 完整 的 Maple 程 序 , 和 典型 问题 的 
数值 计算 数据 , 这 里 的 程序 和 数据 便于 后 续 者 复 现 
计算 和 比较 使 用 . 在 第 3 节 , 对 不 同 参数 的 变化 对 于 
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吸 能 能 力 的 影响 进行 了 细致 的 数值 分 析 . tE SA, 9 
致 计算 分 析 了 吸 能 器 动力 学 的 相 图 在 不 同 参数 变化 
影响 下 的 行为 . 最 后 , 总 结 了 7 点 结论 , 其 中 特别 是 确 
A f Zhang-YuU! 认为 第 开 类 结构 是 “质量 敏感 性 ”的 
科学 论断 
1 折 板 吸 能 器 理论 模型 和 Zhang-Yu 控制 方 
程 图 2 第 I 类 结构 的 四 饮 链 模型 
Fig. 2 Four-hinge model of type II structure 
X 1 ` 5 X " " 
Seque T= x21? 4 Gsin 9 (1) 
Table 1 symbol table 3 
a 系统 在 任 一 时 刻 的 势能 为 = 0， 所以， 系统 
idi a 的 Lagrange( 拉 格 朗 日 ) 是 L = T-V = $mD8 + 
T 系统 在 任 一 时 刻 的 动能 272(m+G)sin? 08^. 由 于 虚 扭 转 66 产生 的 虚 功 为 6W = 
_0 时 刻 系统 动能 Mm 
" de -AM, + M36, 使 用 系统 的 广义 荷载 是 0 = M = 
Uo 初始 输入 动能 
—4(Mı + M2) 
: iR 由 拉 格 朗 日 方程 4 站 - 35 = o 可 以 推导 出 系统 
fo 初始 折 角 BN aab 9 = G 
ò 角速度 运动 的 微分 方程 6 
to t=0 时 刻 角速度 
Mo 准 静 态 完全 塑性 弯 逢 二 L?(m + G) sin cos 002 + Mi+M _ Q) 
70 VEREA JERA L2 [im + (m + G) sin? J 
Mı, M2 (Ez i pc 709 35 2J] 7] B. 
Vo 落 锤 碰撞 初速 度 此 二 阶 微分 方程 的 初始 条 件 可 由 在 冲击 载荷 情况 下 
m 折 板 质量 的 拉 格 朗 日 方程 得 到 , 即 
G TE HEJ EŒ GV, sing 
, 能 量 吸收 系数 Ù = 一 一 一 G) 
试 样 的 宽度 2L |m/3 + (m + G) sin? 6] 
à Eius 中 go 和 多 分别 是 := 0 时 刻 的 折 角 6 和 角速度 6 
3 MED 值 . 这 里 的 公式 的 推导 取 自 Zhang -YuDl, 将 方程 (3) 称 
p 材料 的 密度 


R 


3 HjZhang-YuU! rp F5 Vu ze p 709 (HA 2)— 28 TL 2C 
件 可 以 简化 成 通过 塑性 贸 连 接 的 四 根 刚 性 杆 , 塑性 
贸 处 的 弯 条 Mi 和 M2 可 视 为 主动 力 算 .其 中 L 和 mm 分别 
为 每 根 杆 的 长 度 和 质量 , 所 有 参量 含义 见 表格 1. 假 
设 初 始 速度 为 VW 的 刚体 G 撞 击 样 品 的 顶部 (图 2) 中 的 
点 A, 并 且 在 撞击 后 一 直 沿 点 A 向 下 移动 .相对 应 的 公 
式 如 下 : 


为 Zhang-Yu 方 程 . 

参考 Symonds 等 提出 的 Cowper-Symonds 本 构 关 
系 , Zhang-Yu 中 推导 出 折 板 的 刚 粘 塑 性 模型 和 刚 塑 性 
模型 的 动态 弯 矩 表达 式 如 下 : 


刚 粘 塑性 模型 : 
本 | + (1.94 x 10)*] 0-0 
M= | (4) 
« Mo 当 g =0 
m | 十 (3.88 x 104) 当 g > 0 
M; = (5) 
« Mo 当 g =0 


刚 塑 性 模型 : 
Mo M > 0 
Mı =M = | (6) 
< Mo 140-0 


其 中 Mo = bj2oo/4 是 准 静 态 完 全 塑性 弯 矩 , oo 是 准 静 
态 的 届 服 应 力 . 


2 折 板 吸 能 器 问题 的 Maple 程 序 和 上 典型 计算 
结果 数据 


由 于 Zhang-Yu 方 程 具有 强烈 的 非 线 性 特性 , 其 
线性 化 的 方程 不 能 提供 任何 有 价值 的 物理 信息 , 解 
析 求 解 不 可 能 , 所 以 本 文 将 对 Zhang-Yu 方 程 进 行 数 
HEM. 为 了 便于 结果 验证 和 今后 读者 跟 进 复 现 这 
里 的 研究 , 本 文 特别 使 用 软件 Maple 编 写 了 通用 程序 ， 
数值 计算 使 用 了 Maple 的 4 阶 Runge-Kutta 方 法 的 求解 
器 rkf45, 使 得 求解 这 个 问题 可 以 到 达 很 高 的 精度 , 同 
时 避免 了 读者 编写 程序 的 不 适用 的 问题 . 

从 文献 看 , Zhang-YuD] 文 献 中 给 出 了 部 分 数值 结 
R, 但 没有 给 出 转角 90(?) 的 结果 . m T f fenex 
系统 的 基本 未 知 量 , 有 必要 给 出 详细 的 计算 和 结 
为 了 便于 读者 自己 计算 和 应 用 这 里 的 结果 , 特地 提 
供 刚 粘 塑 性 模型 A 组 的 完整 程序 如 表 2, 其 余 六 组 改 
变相 应 G, Vo 取 值 即 可 . 
相关 系统 参数 与 文献 43 中 所 使 用 的 数值 相同 . 
这 些 是 
样品 的 宽度 , b = 51[mm]; 
试 样 的 厚度 , h = 1.6[mm]; 
样品 的 有 效 长 度 , 2L = 50[mm]; 
单 块 折 板 重量 , m = 0.016[kg]; 
准 静 态 完 全 塑性 弯 矩 , Mo = 9.79[N - m]; 
铁 锤 的 初始 动能 ，Uo = 1220], R34 E T kE 
质量 G 和 碰撞 速度 Vo 七 种 组 合 . 根据 两 类 塑性 理论 ， 
bo = 1.146",2.6",4 三 种 初始 折 角 , 分 为 六 组 进行 分 


表 3 计算 数据 取 自 CALLADINE-ENGLISHOI 
Table 3 Calculated data are taken from 


CALLADINE-ENGLISH!! 


od G Drop height Vo 
(kg) (m) (ms-!) 
A 4.55 2.74 7.33 
B 6.88 1.81 5.96 
C 9.83 1.27 5.00 
D 13.8 0.907 4.22 
E 18.9 0.660 3.60 
F 26.1 0.479 3.07 
G 35.1 0.355 2.64 


为 了 便于 读者 比较 ， 现 给 出 A 组 G = 
4.55[kg], Vo = 7.33[m .ss 情况 下 计算 数据 , 在 20 组 
不 同时 刻 折 角 9 和 6 的 对 应 值 , 详细 数据 见 表格 4. 


3 各 种 参数 变化 对 折 板 吸 能 器 的 吸收 能 量 能 
力 的 影响 
3.1 变形 角度 对 能 量 吸 收 的 影响 


为 了 更 好 对 比 初 始 缺陷 敏感 性 , 在 Zhang-YuD 两 
组 初始 折 角 bp = 1146,80) = 4 的 基础 上 , 增设 一 
组 go = 2.6". 不 同 组 合 的 折 角 6 随时 间 t 的 关系 图 如 下 : 


0 0.0020.0040.0060.0080.0100.012 


A B c D E F G 


图 3 表 3 中 七 组 G,Wo 对 应 数值 下 , 初始 折 角 6 = 1.146" 时 刚 粘 塑性 
Jioc) 图 像 


Fig. 3 0(1) image of rigid-viscoplastic model with initial folding angle 
0o = 1.146? under the corresponding values of seven groups of G and Vo 


in Table 3 


0 0.01 0.02 0.03 
t 


图 4 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 6 = 1.146° 时 刚 塑 性 模 
型 9(7) 图 像 


Fig. 4 0(1) image of rigid-plastic model with initial folding angle 
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表 2 刚 粘 塑性 模型 A 组 程序 示例 


Table 2 Example of group A program for rigid viscoplastic model 


刚 粘 塑性 模型 (程序 ) 
赋值 


L := 25/1000 


m := 0.016 
G := 4.55 


v0 := 7.33 


M0 := 9.79 
本 构 模 型 


conver piecewise (: < 0,0, x 20, M-O- ( 十 ( 


1.94 


3.88 


微分 方程 


1000 | 3] , Heaviside 
M.1 := M. OHeaviside (diff(6(f), t)) + 2.807650451 Heaviside (diff(6(t), r))(diff(9(r), r)) /? 
conver piecewise (: < 0,0, x 20, M-O- ( 十 ( 


l 
m x)? J n Heaviside ) 
M.2 := M_OHeaviside (diff(&(t), t)) + 3.225143455 Heaviside (difr(6(D),D)(difr(6(0D,D)U5 


a := L? - (m + G) - sin(0(®) - cos(6(r)) - diff(G(?), t) - diff(6(r), t) + M1 + M2 
初始 条 件 % = 1.146° 


b := L? - (8 + (m + G) .sin(0D) - sin(6()) 


sys := diff(6(t),t,t) + $ =0 
00) = 146x 


180 > DX8X(0) = 


bes: 
解 方程 


G-v-0-sin( Len) 
4 


2 人 Sm+(m+G) sin( 
方程 


1.14 


180^ y) 
sol : = dsolve (( sys, bcs}, numeric, output —listprocedure) 


im 
Hi, 


122 
gj:1- H 
91 := plots:-odeplot ( sol, [t, 0(t)], 


T:-2.m.LD?-(dift(9(), D)? + 2- L? - (m + G)- (sin(6(0)) - (diff(9(0), D)? 


初始 条 件 g0 = 4? 


t = 0..0.89067682e — 3, color = 
h1 := plots:-odeplot (sol, [t, H], t = 0..0.89067682e — 3, color = red, legend-A) 

f1 := plots:-odeplot( sol, [t, 7], t = 0..0.89067682e — 3, color 
解 方程 


pcs 2 := 6(0) = $4, D(0)(0) = 


Gv_0sin( pza j 
2: 


E 3 m+(m+G) sin( $8) ) 
sol.2 : = dsolve ({ sys, bcs_2}, numeric, output =listprocedure) 
g2 := plots:-odeplot ( sol. 2, [1, 6(1)], 


D 


初始 条 件 6 = 2.6? 


t = 0..0.51093006e - 2, color 


h2 := plots:-odeplot (sol. 2, [t, H], t = 0..0.51093006e — 2, color = blue, legend-A) 
f2 := plots:-odeplot( sol.2, [t, 7], t = 0..0.51093006e — 2, color 
解 方程 


bcs.3 := ((0) = 282, D(9)(0) = 


p^. 


TA 


D 


blue, legend-A) 
G-v- 0-sin( za ) 


22{ 1 m+(m+G) sin( 258) ) 
sol.3 : = dsolve (( sys, bcs_3}, numeric, output =listprocedure) 
g3 := plots:-odeplot ( so1.3, [t, 6(1)], 


t = 0..0.35142103e - 2, color = black, legend = 6) 
h3 := plots:-odeplot (sol.3, [t, H], t = 0..0.35142103e - 2, color = black, legend- A) 


f3 := plots:-odeplot( sol.3, [t,7]], t = 0..0.35142103e — 2, color = black, legend-A) 


red, legend = 6) 


red, legend-A) 


blue, legend = 0) 


Table 4 When G = 4.55[kg], Vo = 7.33[m - s^! ], calculated results data table in group A 


d 4 A 组 G = 4.55[kg], Vo = 7.33[m - s-!] 计 算 结果 数据 表 


0o = 1.146? under the corresponding values of seven groups of G and Vo 


图 5 33H 


in Table 3 


00.002 0.006 0.010 0.014 0.018 
t 


4 


B e D 


E F G 


FP 七 组 G,Vo 对 应 数值 下 , 初始 折 角 6 


型 g(t) 图 


像 


2.6" 时 刚 粘 塑性 


Fig. 5 0(1) image of rigid-viscoplastic model with initial folding angle 


0o = 2.6? under the corresponding values of seven groups of G and Vo 


in Table 3 


图 6 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 


型 6D 图 像 


409 = 2.6? 时 刚 塑 性 模 


Fig. 6 0(1) image of rigid-plastic model with initial folding angle 


0o = 2.6? under the corresponding values of seven groups of G and Vo 


in Table 3 


9 0.5 


0 


刚 粘 塑性 模型 mw = 1.146? 刚 粘 塑性 模型 % = 4 刚 塑性 模型 0 = 1.146? 刚 塑 性 模型 go = 4 
t[ms] 0 O[rad/ s] tims] 0 O[rad/ s] t[ms] 0 O[rad/ s] t[ms] 0 O[rad/ s] 
0 0.0200 1863.2881 0 0.0698 1688.8352 0 0.0200 1863.2881 0 0.0698 1688.8352 
0.0468 0.0734 760.5072 0.2689 0.2583 418.1586 0.1595 0.1401 428.5202 0.7424 0.4265 267.3300 
0.0936 0.1021 505.3475 0.5308 0.3393 289.8303 0.3190 0.1947 280.9162 1.4848 0.5876 180.9920 
0.1404 0.1226 383.7338 0.8067 0.4142 213.2657 0.4785 0.2337 214.1282 2.2272 0.7060 141.7471 
0.1872 0.1387 308.8665 1.0756 0.4658 173.2901 0.6380 0.2644 173.5131 2.9696 0.8017 117.6900 
0.2340 0.1519 256.4686 1.3445 0.5084 145.0008 0.7975 0.2897 145.1565 3.7120 0.8824 100.6927 
0.2808 0.1629 216.8262 1.6134 0.5444 123.4842 0.9570 0.3110 123.6722 4.4544 0.9521 87.6050 
0.3276 0.1723 185.2092 1.8823 0.5752 106.1670 1.1165 0.3293 106.4785 5.1968 1.0131 76.9165 
0.3744 0.1803 159.0061 2.1512 0.6017 91.6778 1.2760 0.3452 92.1613 5.9392 1.0667 67.8005 
0.4212 0.1872 136.6418 2.4201 0.6246 79.1916 1.4355 0.3588 79.8710 6.6816 1.1140 59.7620 
" 0.4680 0.1931 117.1026 2.6890 0.6444 68.1783 1.5950 0.3707 69.0607 7.4240 1.1557 52.4841 
0.5148 0.1982 99.7020 2.9579 0.6614 58.2798 1.7545 0.3809 59.3591 8.1664 1.1921 45.1540 
0.5616 0.2025 83.9569 3.2268 0.6758 49.2445 1.9140 0.3897 50.5020 8.9088 1.2237 39.4231 
- 0.6084 0.2061 69.5164 3.4957 0.6880 40.8904 2.0735 0.3971 42.2945 9.6512 1.2507 33.3847 
0.6552 0.2090 56.1206 3.7646 0.6979 33.0834 2.2330 0.4032 34.5872 10.3936 1.2733 271.5605 
0.7020 0.2114 43.5761 4.0335 0.7058 25.7243 2.3925 0.4081 27.2619 11.1360 1.2917 21.8917 
0.7488 0.2130 32.7292 4.3024 0.7118 18.7425 2.5520 0.4119 20.2218 11.8784 1.3058 16.3333 
0.7956 0.2143 20.5345 4.5713 0.7159 12.0960 2.7115 0.4146 13.3846 12.6208 1.3159 10.8502 
0.8424 0.2150 9.9420 4.8402 0.7183 5.7864 2.8710 0.4162 6.6780 13.3632 1.3220 5.4138 
0.8906 0.2153 0.0125 5.1090 0.7191 0.0052 3.0313 0.4167 0.0020 14.1050 1.3240 0.0042 


孙 博 华 等 :Caladine-English 折 板 吸 能 器 的 能 量 吸 收 能 力 分 析 


图 7 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 00 
型 9(7) 图 像 


4 时 刚 粘 塑性 模 


Fig. 7 0(1) image of rigid-viscoplastic model with initial folding angle 


0o = 4? under the corresponding values of seven groups of G and Vo in 


Table 3 


图 8 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 60 
型 g(z) 图 像 


Fig. 8 0(1) image of rigid-plastic model with initial folding angle 


4 时 刚 塑性 模 


00 = 4? under the corresponding values of seven groups of G and Vo in 


Table 3 


图 例 中 A、B、C、D、E、F、G 分 别 对 应 于 表 2 中 G, Vo 取 
值 情况 . 以 图 3 中 A 组 蓝 线 为 例 进 行 说 明 : 图 中 蓝 线 表 
示 G=4.55 [kg], Vo = 733[m - s-1]( 见 表 3 中 A 组 数据 )， 
采用 刚 粘 塑性 理论 模型 , 在 初始 折 角 6 = 1.146* 时 计 
算 所 得 到 的 折 角 69 与 时 间 t 的 关系 图 . 从 以 上 图 形 可 以 
得 到 以 下 结论 : 

(025 7E fi i A 2 BEUS-122[J]79 XE fü EF, — BG EE 
度 m 逐 渐 减 小 ( 落 锤 重量 G 相 应 增 大 )， 最 终 折 角 9 增 
K, 即 碰撞 损失 的 能 量 减 少 . 这 与 文献 42303133 中 所 
得 到 的 第 开 类 结构 “速度 /质量 敏感 性 ”的 结论 是 完 
一 致 的 . 
CO) 随 着 初始 折 角 6 的 增 大 , 无 论 是 刚 粘 塑性 模型 还 
是 刚 塑性 模型 , A-G 组 中 相 邻 组 之 间 最 终 折 角 6 的 差 
值 明 显 减 小 , 质量 敏感 性 减弱 ,而 且 结构 变形 的 速 
度 更 加 缓慢 . 这 与 Zhang-YuD 中 强调 的 “初始 几何 缺 
陷 ” 对 “质量 敏感 性 ”的 显著 影响 是 吻合 的 . 

(3) 当 go = 1.146" 时 , 尽管 速度 VW 不断 改变 , 其 快速 变 
形 阶段 变形 速率 保持 一 致 ; 但 是 随 着 初始 折 角 的 
增 大 , 不 同 组 间 变 形 速 度 随 着 Vo 改变 明显 分 层 , 初 
始 速 度 Vo 越 大 (G 相 应 越 小 ), 变形 速度 越 快 . 余 同 希 
等 09 给 出 了 解释 : 当初 始 折 角 较 小 时 , 横向 惯性 对 折 
板 吸 能 器 的 动力 行为 起 支配 作用 . 
(4) 与 刚 塑性 模型 相 比 , 刚 粘 塑性 模型 
降低 , 且 其 变形 过 程 更 加 迅速 . 这 可 能 
流动 与 应 变 率 效 应 导致 的 . 


最 终 转角 大 幅 
是 考虑 了 塑性 


32 ”质量 比 和 初始 速度 对 能 量 吸收 的 影响 


在 系统 响应 的 第 二 相 , 刚性 杆 稳定 地 绕 着 四 个 
铵 转动 ,其 最 终 转 角 与 系统 第 一 相 碰撞 后 剩余 的 
能 量 70 成 正比 .Zhang-YuDI 推 导出 公式 (89) 并 指出 ， 
Uo = GV3/2 不 变 时, 初速 度 Vo 的 改变 是 由 质量 G 的 
改变 导致 的 . 为 了 进一步 看 出 初始 速度 V6。 和 To 之 间 
的 关系 . 由 公式 (8) 绘 出 了 7o/Vo-G/m 和 7o/Vo-W 的 
图 像 . 图 中 实 线 、 虚 线 、 点 划 线 分 别 代表 初始 折 角 go 
等 于 1.146"、2.6*、4° 的 情况 . 


1000 3000 5000 7000 10000 
GS 
m 
[I— 16 ——26-:-4 


图 9 7o/Vo 与 G/m 的 关系 图 


Fig. 9 The diagram of 70/Vo and G/m 
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图 10. To/U6o 与 初始 速度 Vo 的 关系 图 


Fig. 10 The diagram of 70/Vo and the initial velocity Vo 


从 上 图 可 以 看 出 ， 

(在 碰撞 的 瞬时 有 部 分 能 量 损 失 , 质量 比 m/G 越 大 ， 

初始 速度 Vo 越 大 , 能 量 损失 的 越 多 ， 

(CO) 初 始 缺陷 对 质量 敏感 性 影响 十 分 显著 . 随 着 初始 

折 角 6 的 增 大 , 质量 敏感 性 极 大 削弱 , 而 且 随 着 质量 

比 的 增 大 (Vo 相应 的 增 大 ), 削弱 的 程度 越 明显 . 
为 了 进一步 了 解 系统 响应 的 第 二 相 的 能 量 耗 散 

情况 , 将 求解 出 的 数值 解 代 入 方程 (D) 绘 出 了 如 下 情 

况 的 7/Vo — ! 的 关系 图 . 以 下 图 形 根据 三 种 不 同 初始 

折 角 值 和 两 种 塑性 理论 分 为 六 组 : 


E 
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DE 
0.002 0.006 0.010 


图 11 表 3 中 七 组 G,Wo 对 应 数值 下 , 初始 折 角 60 = 1.146° 时 刚 烙 塑性 
EU TOR 


Fig. 11 Ln (t) image of rigid viscoplastic model with initial folding 
angle 69 = 1.146? under the corresponding values of seven groups of G 


and Vo in Table 3 


图 12 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 00 = 1.146°? 时 刚 塑 性 模 
0 T OE 


Fig. 12 5 (t) image of rigid plastic model with initial folding angle 
0o = 1.146? under the corresponding values of seven groups of G and Vo 


in Table 3 


0 
0 0.006 0.010 0.014 0.018 
t 
4 B E D E F G 


图 13 表 3 中 七 组 G,W 对 应 数值 下 , 初始 折 角 00 = 2.6^ 时 刚 粘 塑性 模 
H ORNS 


Fig. 13 Ln (t) image of rigid viscoplastic model with initial folding 


angle 00 = 2.6? under the corresponding values of seven groups of G 


and Vo in Table 3 


图 14 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 b = 2.6" 时 刚 塑 性 模 
m TO 图 像 


Fig. 14 5,0 image of rigid plastic model with initial folding angle 
0o = 2.6? under the corresponding values of seven groups of G and Vo 


in Table 3 


0.6 


0.005 0.010 0.015 0.020 
t 
A B c D E F G 


图 15 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 bo = 4° 时 刚 粘 塑 性 模 
H TO 图 像 


Fig. 15 i0 image of rigid viscoplastic model with initial folding 
angle 00 = 4° under the corresponding values of seven groups of G and 


Vo in Table 3 


0 001 002 003 004 


图 16 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 6 = 4 时 刚 塑 性 模 
MEO 图像 


Fig. 16 5,0 image of rigid viscoplastic model with initial folding 
angle 00 = 4° under the corresponding values of seven groups of G and 


Vo in Table 3 


图 例 中 A、B、C、D、E、F、G 分 别 对 应 于 表 3 中 G, Vo W 
值 情 况 . 以 图 11 中 A 组 蓝 线 为 例 进行 说 明 : 图 中 蓝 线 
表示 G = 4.55[kg], Vo = 7.33[m:s-!]( 见 表 3 中 A 组 数据 )， 
采用 刚 烙 塑性 理论 模型 , 在 初始 折 角 go = 1.146°* 时 计 
算 所 得 到 的 能 量 比 T/U6 与 时 间 t 的 关系 图 . 从 上 列 图 
形 中 可 以 看 出 : 
(1t 等 于 0 的 时 刻 , To 的 值 与 图 9 保持 一 致 , 且 刚 粘 塑性 
模型 与 刚 塑 性 模型 没有 差异 . 

(2) 质 量 比 m/G 越 大 (初始 速度 Wo 越 大 ), 系统 响应 的 第 
二 相 耗 能 越 快 , 而 初始 折 角 go 减缓 这 种 趋势 . 

(3) 刚 粘 塑 性 模型 与 刚 塑 性 模型 耗 能 趋势 总 体 保 持 一 
致 , 后 者 耗 能 时 间 大 约 是 前 者 的 两 到 三 倍 . 
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3.3. 折 板 吸 能 器 的 吸收 能 量 能 力 分 析 


为 了 定量 的 预测 结构 的 吸 能 能 力 , 引入 一 个 能 
量 吸收 系数 
HO qe 
w=- cam (7) 


0+ 
0.002 0.004 0.006 0.008 0.010 0.012 
t 
4 B Ẹ D E F G 


图 17 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 go = 1.146" 时 刚 粘 塑性 
模型 能 量 吸收 系数 7D 图 像 


Fig. 17 r(t) image of rigid viscoplastic model with initial folding angle 


0o = 1.146? under the corresponding values of seven groups of G and Vo 


in Table 3 


0.01 0.02 0.03 


E 
a 
[ 
G 
^ 
a 


图 18 表 3 中 七 组 G,Vo 对 应 数值 下 , WRIT = 1.146° 时 刚 塑 性 模 
型 能 量 吸收 系数 7(D 图 像 


Fig. 18 7(t) image of rigid plastic model with initial folding angle 
0o = 1.146? under the corresponding values of seven groups of G and Vo 


in Table 3 


t 


A B € D E F G 


图 19 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 6 = 2.6? 时 刚 粘 塑性 模 
型 能 量 吸收 系 数 7(D 图 像 


Fig. 19 ņ(t) image of rigid viscoplastic model with initial folding angle 
0o = 2.6? under the corresponding values of seven groups of G and Vo 


in Table 3 


DS 


20 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 90 = 2.6^ 时 刚 塑 性 模型 
能 量 吸 收 系数 7(D 图 像 


Fig. 20 7(D image of rigid plastic model with initial folding angle 
0o = 2.6? under the corresponding values of seven groups of G and Vo 


in Table 3 


0+ 

0 0.005 0.010 0.015 0.020 
t 

A B [3 D E F G 


图 21 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 6 = 4? 时 刚 粘 塑性 模型 
能 量 吸收 系数 7(D 图 像 


Fig. 21 7(t) image of rigid viscoplastic model with initial folding angle 
0o = 4? under the corresponding values of seven groups of G and Vo in 


Table 3 


IR] 


22 表 3 中 七 组 G,Vo 对 应 数值 下 , 初始 折 角 go = 4° NABER fe 
量 吸收 系数 (7) 图 像 


Fig. 22 n(t) image of rigid plastic model with initial folding angle 
0o = 4? under the corresponding values of seven groups of G and Vo in 


Table 3 


图 例 中 A、B、C、D、E、F、G 分 别 对 应 于 表 3 中 G, Vo 取 
值 情况 . 以 图 17 中 A 组 蓝 线 为 例 进行 说 明 : 图 中 蓝 线 
表示 G = A 55[kg], Vo = 7.33[m:s-!]( 见 表 3 中 人 A 组 数据 )， 
采用 刚 粘 塑性 理论 模型 , 在 初始 折 角 6 = 1.146 B] TF 
算 所 得 到 的 能 量 吸收 系数 7 与 时 间 t 的 关系 图 . 

从 上 列 图 形 中 可 以 看 出 : 随 着 质量 G 的 增 大 ( 初 
速度 Ww 相应 减 小 ), 能 量 吸收 系数 7 相应 减 小 ; 初始 折 
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角 % 对 能 量 吸收 系数 ”影响 显著 , % 越 小 , 对 7 影响 越 
显著 , 这 与 之 前 得 到 的 结论 是 一 致 的 . 


对 于 时 刻 ; - 0, 联 立 方程 (1) 和 (3)， 
Zhang 和 Yu 中 得 到 
To m 1 
于 | 未 = 1 
Uo | 和 4 3 n) < " 


所 以 ,t= 0 时 刻 的 能 量 吸收 系数 是 xx(0) = 1- Dj 


m 1 
G ( F s) 


Le 81s 


1(0) = (9) 


3sin? i] 


RPU = S ict M HO GE, To 是 碰撞 后 瞬 
间 系 统 动能 , 即 公式 (9) PTE = 0 时 的 值 . 


图 23 公式 (9)mo - 6 - 各 三 维 图 


Fig. 23 o — 6o — c three-dimensional diagram of formula (9) 


71(0) = 1 一 To/Uo 表 示 在 初始 时 刻 , 系统 有 能 量 损失 ， 
初始 角度 0 的 影响 很 小 , 质量 比 m/G 的 影响 比较 显 
著 , 既 杆 的 重量 m 越 大 , 吸收 刚体 质量 G 冲 击 的 能 量 
能 力 越 强 . 


0 
0.006,0006 0.0003 


图 24 刚 粘 塑性 A 组 bo = 1.146" 时 7(D - 00) — t — 2 i £s End 


Fig. 24 Three-dimensional graph of 77(D — 6(t) — t for rigid 


viscoplasticity group a with 69 = 1.146? 


4 转速 动力 学 相 图 分 析 
41 初始 缺陷 变化 的 动力 学 相 图 分 析 
为 了 定量 分 析 初 始 缺陷 对 于 吸 能 结构 的 影响 ， 


以 A 组 数据 为 例 , 下 面 给 出 了 刚 粘 塑性 和 刚 塑 性 模型 
下 不 同 初 始 折 角 的 0 — 9 POSER ES. 
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25 刚 粘 塑 性 模型 A 组 : 不 同 初始 折 角 go = 10/2 F0- d fcx 
相 图 , 其 中 i=1,2...10 


Fig. 25 Rigid viscoplastic model group A: 0 — e rotational speed phase 


DS 


diagram with different initial folding angles 0o = i?/2, where i-1,2...10 
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图 26 刚 塑 性 模型 A 组 : 不 同 初始 折 角 60 = zz/2 下 0- EFRR, H 
中 i=12..10 


Fig. 26 Rigid plastic model group A: 0 — 4r 2 rotational speed phase 


diagram with different initial folding angles 00 = i? /2, where i-1,2...10 


其 中 i=1,2...10, 相对 应 的 曲线 初始 折 角 为 p = P2. 
以 图 25 中 蓝 线 i=1 为 例 进行 说 明 : 图 中 蓝 线 表示 G = 
4.55[kg], Vo = 7.33[m - s-!]( 见 表 3 中 A 组 数据 ), 采用 刚 
粘 塑 性 理论 模型 , 在 初始 折 角 6 = 0.5° 时 计算 所 得 到 
的 转角 6 与 角速度 风 的 转速 相 图 . 

从 以 上 相 图 可 以 看 出 : (GD 刚 粘 塑性 模型 和 刚 塑 
性 模型 计算 出 的 在 相 平面 中 的 运动 轨迹 大 致 相同 ; 
(2) Eo < 2 时 , 随 着 初始 折 角 的 增 大 , 初始 角速度 
和 最 终 折 角 值 也 增 大 (碰撞 损失 能 量 减少 ), 这 与 之 前 
得 到 的 结论 一 致 ; (3)go = 0.5° 和 go = 1° 时 , 折 板 吸 能 
器 运动 轨迹 与 其 他 各 组 存在 明显 差异 , 此 时 初始 角 
速度 大 幅 降低 (碰撞 损失 能 量 增 大 ), 折 板 吸 能 器 吸 能 
能 力 优异 , 这 也 能 更 好 看 出 初始 缺陷 对 折 板 吸 能 
动态 行为 的 显著 影响 . 


42 ”初始 碰撞 速度 变化 的 动力 学 相 图 分 析 


为 了 定量 分 析 初 始 碰撞 速度 Wo 对 于 吸 能 结构 的 
影响 , 在 保证 输入 动能 Do = 122[J] 不 变 的 情况 下 ， 
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下 面 给 出 了 刚 粘 塑性 和 刚 塑 性 模型 下 初始 折 角 分 别 
为 1.146"、4* 时 的 不 同 初始 碰撞 速度 9 - 内 转速 相 图 . 
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P. 1000 
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27 刚 粘 塑 性 模型 组 : eg = 1.1469 TE, 不 同 初 始 碰撞 速 
度 W = ipn: s] Fe- HPEH R, 其 中 i=1,2...10 


Fig. 27 Rigid viscoplastic model: when 6, = 1.146?, 0 — e rotational 


DS 


speed phase diagram with different impact velocities Vo = i[m - s71], 


where i-1,2...10 
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图 28 刚 塑 性 模型 组 : 6o = 1.146? 时 , 不 同 初始 碰撞 速 
度 Vo = ilm- s-!] 下 0- 邮 转速 相 图 , 其 中 i=1,2...10 


i 


Fig. 28 Rigid plastic model group: when 6o = 1.146^, the phase 
diagram of 0 - 7r ? under different initial collision speeds Vo = i[m - s EUR 


where i-1,2...10 


图 29 刚 粘 塑性 模型 组 : go = 4^ TP, 不 同 初始 碰撞 速 
度 W = ilm- s-!] 下 9 一朗 转速 相 图 , 其 中 i=1,2...10 


Fig. 29 Rigid viscoplastic model: when 60 = 4°, 0 — do rotational speed 
phase diagram with different impact velocities Vo = i[m - s^! ], where 


iz1,2...10 


图 30 刚 塑 性 模型 组 : bo = 4 时 , 不 同 初始 碰撞 速度 Wo = im - s7!] 
Fo- 多 转速 相 图, 其 中 i=1,2...10 


Fig. 30 Rigid plastic model: when 60 = 4°, 0 — 多 rotational speed 
phase diagram with different impact velocities Vo = i[m - s71], where 


iz1,2...10 


| 


其 中 i=1,2...10, 相对 应 的 曲线 初始 碰撞 速度 为 V。= 

1. 以 图 27 中 蓝 线 i=1 为 例 进行 说 明 : 图 中 蓝 线 
表示 落 锤 初始 碰撞 速度 Vo = 4[m:s-!], 采用 刚 粘 塑性 
里 论 模 型 , 在 初始 折 角 go = 1.146* 时 计算 所 得 到 的 转 
fho 与 角速度 蝶 的 转速 相 图 . 

从 以 上 相 图 可 以 看 出 : (1) 当 初始 折 角 很 小 时 , 相 
平面 折 板 吸 能 器 运动 轨迹 较为 混乱 , 随 着 初始 折 角 
增 大 , 相 图 分 布 趋 于 规律 且 运 动 趋势 类 似 ; QT A 
撞 速 度 很 小 时 , 速度 的 改变 对 折 板 吸 能 器 吸 能 能 
影响 显著 , 且 随 着 初始 碰撞 速度 的 增 大 折 板 吸 能 器 
吸 能 能 力 减弱 . 


n 


| 


43 落 锤 质量 G 变 化 的 动力 学 相 图 分 析 


为 了 定量 分 析 落 锤 质 量 对 于 吸 能 结构 的 影响 
在 保证 输入 动能 VU。 = ”122[J 四 不 变 的 情况 下 ， " 
面 给 出 了 刚 粘 塑性 和 刚 塑 性 模型 下 初始 折 角 分 别 
为 1.146"、4" 时 不 同 落 锤 质量 的 0 — 几 转 速 相 图 . 
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31 刚 粘 塑性 模型 组 : go = 1.146? E, 不 同 落 锤 质量 G = Ai[kg] 
Fo- ZHR, 其 中 i=1,2...10 


Fig. 31 Rigid viscoplastic model: when 00 = 1.146°,0— de rotational 
speed phase diagram with different drop weights G = Ai[kg], where 


iz1,2...10 
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图 32 刚 塑性 模型 组 : 6o = 1.146? 时 , 不 同 落 锤 质 
ÆG = Ai[kg] F0 — 内 转速 相 图 , 其 中 i=1,2...10 


Fig. 32 Rigid plastic model: when 00 = 1.146°, 0 — de rotational speed 


phase diagram with different drop weights G = 4i[kg], where i-1,2...10 
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33 刚 粘 塑性 模型 组 : go = 4 时 , 不 同 落 锤 质量 G = 4i[ke] 
Fo- d eSERBES, 其 中 i=1,2...10 


Fig. 33 Rigid viscoplastic model: when 60 = 4°, 0 — do rotational speed 


phase diagram with different drop weights G = 4i[kg], where i=1,2...10 
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图 34 刚 塑性 模型 组 : 90 = 4° 时 , PEHE EG = 4ilkg] F0 — £ 
速 相 图 , 其 中 i=1,2...10 


Fig. 34 Rigid plastic model: when 60 = 4?,0 — e rotational speed 


phase diagram with different drop weights G = 4i[kg], where i=1,2...10 


其 中 i=12...10, 相 对 应 的 曲线 落 锤 质量 为 G = Ai [kg]. 
以 图 31 中 蓝 线 i = 1 为 例 进行 说 明 : 图 中 蓝 线 表示 落 
EEG = A[kg] , 采用 刚 粘 塑 性 理论 模型 , 在 初始 折 
fabo = 1.146" 时 计算 所 得 到 的 转角 6 与 角速度 几 的 转 
速 相 图 . 

从 以 上 相 图 可 以 看 出 : GD) 刚 粘 塑性 模型 和 刚 塑 
性 模型 计算 出 的 在 相 平面 中 的 运动 轨迹 大 致 相同 ; 
Q) 当 落 锤 质 量变 化 时 , 初始 折 角 很 小 的 情况 下 , 折 
板 吸 能 器 运动 轨迹 十 分 相似 , MWER, 即 此 时 折 板 
吸 能 器 吸 能 能 力 变 化 很 小 . 随 着 初始 折 角 增 大 , 相 图 


旺 规律 分 布 , 此 时 随 着 落 锤 质量 增 大 折 板 吸 能 器 吸 
能 能 力 显著 增强 . 


44 折 板 质量 m 变 化 的 动力 学 相 图 分 析 


为 了 定量 分 析 折 板 质量 m 对 于 吸 能 结构 的 影 
响 , 在 保证 输入 动能 Do = 122[J 力 不 变 的 情况 下 ,下 
面 给 出 了 刚 粘 塑性 和 刚 塑性 模型 下 初始 折 角 分 别 


为 1.146*、4° 时 不 同 折 板 质量 的 9 — 学 转速 相 图 . 
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刚 粘 塑性 模型 A 组 : g = 1.1469 FE, 不 同 折 板 质 
tm = 0.01. ilke] FO- 22 转速 相 图 , H xe 


Fig. 35 Rigid viscoplastic model: when 69 = 1.146^, 


X 


lm 


= ? rotational 
speed phase diagram with different drop weights m = 0.01 - i[kg], where 
iz1,2...10 


图 36 刚 塑 性 模型 A 组 : 09 = 1.146° 时 , 不 同 折 板 质 
量 m = 0.01.i[kg] 下 9- 9? 转速 相 图 , 其 中 i=1,2...10 


Fig. 36 Rigid plastic model: when 00 = 1.146°, 


- rd rotational speed 
phase diagram with different drop weights m = 0.01 - i[kg], where 


iz1,2...10 


DS 


37 刚 粘 塑性 模型 A 2H: 99 = A* ENT, 不 同 折 板 质量 mm = 0.01 - i[kg] 
Fo- 29 转速 相 图 , 其 中 i=1,2...10 


Fig. 37 Rigid viscoplastic model: when 00 = 4°, 6 — d rotational speed 
phase diagram with different drop weights m = 0.01 - i[kg], where 


iz1,2...10 
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= i=2 i4 m 
i=6 i=7 ies i=9 i=10 


图 38 刚 塑性 模型 A 组 : b = 4 时, 不 同 折 板 质量 m = 0.01 - ifkg] 
Fo- f cit E, 其 中 i=1.2.…10 


Fig. 38 Rigid plastic model: when 00 = 4°, 0 — e rotational speed 


phase diagram with different drop weights m = 0.01 - i[kg], where 


iz1,2...10 


其 中 i=1,2...10, 相 对 应 的 曲线 折 板 质量 为 m = 0.01 - 
ifkg]， 以 图 35 中 蓝 线 i = 1 为 例 进行 说 明 : 图 中 蓝 线 
表示 折 板 质量 m = 0.01[kg] , 采用 刚 粘 塑性 理论 模型 ， 
在 初始 折 角 go = 1.146* 时 计算 所 得 到 的 转角 9 与 角 速 
度 多 的 转速 相 图 . 

从 以 上 相 图 可 以 看 出 : 刚 粘 塑 性 模型 和 刚 塑 性 
模型 计算 出 的 在 相 平 面 中 的 运动 轨迹 大 致 相同 ; 值 
得 注意 的 是 , 当 m < 0.03[kg] 时 , 随 着 折 板 质量 的 增加 ， 
撞 瞬 间 损 失 能 量 迅速 增加 , 折 板 吸 能 器 运动 轨迹 
化 很 大 . 随 着 折 板 质量 继续 增 大 , 相 图 变化 不 明显 ， 
对 折 板 吸 能 器 吸 能 影响 减弱 . 


45 折 板 与 落 锤 质量 比 myG 的 动力 学 相 图 分 析 


为 了 定量 分 析 折 板 与 落 锤 质量 比 m/G 对 于 吸 能 
结构 的 影响 ,在 保证 输入 动能 Uo = 122[J] 不 变 的 情况 
下 , 下 面 给 出 了 刚 粘 塑性 和 刚 塑 性 模型 下 初始 折 角 
分 别 为 1.146"、4* 时 不 同 质量 比 的 0 — 学 转 速 相 图 . 


DS 


39 刚 粘 塑 性 模型 A 组 : 6o = 1.1462 EF, 不 同 质量 
比 m/G = y Fe - 22 转 速 相 图 , 其 中 i=1,2...10 


Fig. 39 Rigid viscoplastic model: when 00 = 1.146°,0— de rotational 
speed phase diagram with different mass ratios m/G — xb where 


iz1,2...10 


图 40 du pape" % = = 1.146° 时 , 不 同 质量 
比 m/G= qd ES, 其 中 i=1,2...10 


Fig. 40 Rigid plastic model: when 00 = 1.146?, 0 — d rotational speed 


phase diagram with different mass ratios m/G — 3b where i-1,2...10 
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图 41 刚 粘 塑性 模型 A 2H: 6 = 4° 时 , 不 同 质量 比 m = 0.01 - i[kg] 
Fo- 4o 转速 相 图 , 其 中 i=1,2...10 


Fig. 41 Rigid viscoplastic model: when 69 = 4°, 0 — 2 rotational speed 


phase diagram with different mass ratios m/G = sin where i-1,2...10 


图 42 刚 塑性 模型 A 组 : 90 = 4? 时 , 不 同 质量 比 m/G = xd; Fo - de 
速 相 图 , 其 中 i=1,2...10 


Fig. 42 Rigid plastic model: when 00 = 4°,0— rd rotational speed 


phase diagram with different mass ratios m/G = b where i-1,2...10 


其 中 i=1,2.….10, 相 对 应 的 曲线 质量 比 为 m/G = agg. 以 
图 39 中 蓝 线 ; = 1 为 例 进行 说 明 : 图 中 蓝 线 表示 质 
量 比 m/G = qs. 采用 刚 粘 塑性 理论 模型 , 在 初始 折 
角 6 = 1.146" 时 计算 所 得 到 的 转角 6 cs fe RE P RO 
速 相 图 . 
从 以 上 相 图 可 以 看 出 : (GD) 刚 粘 塑性 模型 和 刚 塑 
性 模型 计算 出 的 在 相 平面 中 的 运动 轨迹 大 致 相同 , 特 
别 是 初始 折 角 较 大 时 近乎 重 车 . 但是, 当初 始 折 角 很 
小 时 , 随 着 质量 比 的 变化 相 图 明显 分 层 , 即 质量 比 对 
折 板 吸 能 器 吸 能 能 力 影响 显著 ; (2) 当 质量 比 变化 时 ， 
初始 折 角 对 相 图 变化 的 影响 相 较 于 其 他 几 组 更 小 . 
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5 结论 致谢 : 余 同 希 教授 是 国际 结构 弹 塑性 动力 学 等 

PUN 领域 的 著名 学 者 , 本 文 是 基于 余 先生 的 一 项 研究 成 

本 文采 用 工程 计算 软件 Maple 数 值 计算 的 方法 ， 果 四 作者 希望 以 此 文 向 余 先 生 致敬 .同时 作者 还 想 
在 输入 动能 恒定 时 , 对 7 组 不 同 质量 G( 初 始 速度 ) 的 落 


Te 感谢 研究 生 戴 远 帆 , 郭 晓 琳 , 吴 凡 对 论文 中 图 3- 图 8 的 
锤 进行 了 折 板 的 动态 行为 研究 , 重点 分 析 了 不 同 塑性 。 整理 和 提供 的 一 此 宝贵 建议 | 
理论 、 不 同 初始 折 角 、 不 同 质量 G( 初 始 速度 ) 对 折 板 
结构 动态 行为 的 影响 , 得 到 如 下 主要 结论 : a x dx 
(1)* 速 度 敏感 性 ”与 碰撞 时 的 瞬时 能 量 损失 密切 


> 


HX, 随 着 初始 速度 Vo 增 大 , 能 量 损 失 增 大 , 结构 最 
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